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ABSTRACT 

Observations of the CMB can be contaminated by diffuse foreground emission from sources 
such as Galactic dust and synchrotron radiation. In these cases, the morphology of the contam- 
inating source is known from observations at different frequencies, but not its amplitude at the 
frequency of interest for the CMB. We develop a technique for accounting for the effects of such 
emission in this case, and for simultaneously estimating the foreground amplitude in the CMB 
observations. We apply the technique to CMB data from the MAXIMA-1 experiment, using 
maps of Galactic dust emission from combinations of IRAS and DIRBE observations, as well as 
compilations of Galactic synchrotron emission observations. The spectrum of the dust emission 
over the 150-450 GHz observed by MAXIMA is consistent with preferred models but the effect 
on CMB power spectrum observations is negligible. 

Subject headings: cosmic microwave background - cosmology: observations 

1. Introduction 

Measurements of the Cosmic Microwave Background (CMB) have begun to fulfill their promise to 
image the Universe at the epoch of the decoupling of photons from baryons, thereby measuring cosmological 
parameters to new levels of precision. In recent years, the balloon-borne bolometer experiments MAXIMA-1 
(Hanany et al. 2000; Lee et al. 2001; Rabii et al. 2003), BOOMERANG (de Bernardis et al. 2000; Netterfield 
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ct al. 2002) and Archeops (Bcnoit et al. 2002), as well as the ground-based interferometers, CBI (Padin 
ct al. 2001; Mason et al. 2002), DASI (Halverson et al. 2002) and VSA (Grainge et al. 2002) have measured 
the CMB power spectrum down to angular scales of 10 arcminutes or better; most recently, WMAP has 
observed the full CMB sky over a factor of four in frequency and with a beam of 13 arcmin FWHM. The 
CMB photons are produced at the last scattering surface, when the opaque, charged plasma of electrons 
and nuclei becomes a transparent gas of hydrogen and helium, at a temperature of about 1 eV, the epoch 
of "Recombination." The measured signals thus map the primordial distribution of matter in the Universe, 
and allow measurements of the cosmological parameters and an understanding of the mechanism for laying 
down the initial fluctuations. 

This promise has, however, always been tempered by the possibility of contamination from astrophysical 
sources of microwave emission or absorption, even away from the Galactic plane. Happily, the cosmological 
signal has proven to be dominant over a wide range of frequencies and angular scales. Nonetheless, we can 
hope to isolate these foreground contributions: the CMB itself has the well-understood shape of a black-body 
spectrum at a given temperature. Other contributions, Galactic or extragalactic, will in all but the most 
pathological cases have a different spectral dependence. While the peak of the 2.73K CMB intensity lies at 
about 90 GHz, other contributions dominate at other frequencies (e.g., Tegmark et al. 2000). Emission from 
dust (and dusty external galaxies) is expected to dominate the CMB at higher frequencies; synchrotron and 
free-free emission (and galaxies radiating by these methods) dominate at lower frequencies. 

Very often, then, we will have independent measurements of the contribution from a given foreground 
component. However, because those measurements are taken at a different frequency than those of the CMB 
measurements, we must extrapolate down to the frequencies of interest. As we inevitably lack a perfect 
understanding of the foreground emission mechanism, this extrapolation will be imprecise. 

In this work, we present a framework for dealing with such imprecise extrapolation, by leaving the spec- 
tral dependence free to vary, but fixing the morphology to be determined by the external measurements. The 
method allows us, on the one hand, to determine the global spectral behavior of the foreground component 
and "marginalize over" (in Bayesian parlance) the CMB signal, or, on the other hand, to determine the CMB 
signal while marginalizing over the foreground contribution. Similar methods have been proposed in the past 
with more ad hoc derivations (e.g., Dodelson & Stebbins 1994; Tegmark & Efstathiou 1996; Dodelson 1997; 
Tegmark 1998; de Oliveira-Costa & Tegmark 1999, and references therein). Here, we show how the same 
formalism can deal with several different foreground problems: estimating the amplitude of the foreground 
emission; estimating the CMB map after accounting for the foregrounds; and finally estimating the CMB 
power spectrum in the presence of such contamination. 

In this paper we concentrate on the contribution of Galactic dust emission to the emission observed by 
the MAXIMA-1 experiment at 150-410 GHz. As our foreground template, we use the recent seminal work of 
Finkbeiner et al. (1999), who combined data from the IRAS satellite with that from the DIRBE and FIRAS 
instruments on the COBE satellite to extrapolate the spectrum of Galactic dust pixel-by-pixel to our CMB 
frequencies. These maps were also used in the analysis of the dust signal present in BOOMERANG data 
(Masi et al. 2001). 

Finally, we note that the recent work of the WMAP team takes a somewhat different approach, using a 
maximum-entropy method to estimate foreground emission (Bennett et al. 2003). In that work, they too use 
the Finkbeiner et al. (1999) maps as a "prior" for the dust emission. However, the effective sensitivity of the 
WMAP dust reconstruction is comparable to the CMB sensitivity. Thus, in regions of very low dust contrast 
such as that observed by MAXIMA-1, the WMAP prediction of dust emission is insufficiently sensitive, and 
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indeed imperfectly correlated with the input prior maps. 

2. Methods 

In this paper, we are concerned with the possible contribution from the aforementioned sources of 
Galactic foreground emission, in particular contamination by dust emission. We will further assume that 
we have a robust measurement of the foreground at some frequency other than that with which we observe 
the CMB. For example, dust is dominant at high frequencies, / > 300 GHz, but believed to be only a small 
contaminant at the 50-200 GHz at which the CMB dominates high-Galactic latitude emission. 

Thus, we will assume that the spatial pattern of these sources is known, and described by some template, 
fi p , where p numbers the observed pixels in our CMB maps, and i labels the various foreground contributions 
(e.g., dust, synchrotron, etc.). However, we will leave the amplitude of the individual foreground contributions 
to be determined. That is, fi p is the foreground map, at its observed frequency for example. This will be 
multiplied by an unknown amplitude Pi, to be determined as the result of our analysis. 

For the purposes of this paper, then, the intensity in a pixel of our CMB map is given by 

d p = X! & f*P + s p + n P W 

i 

where d p is the data at pixel p; s p is the CMB signal at that pixel; n p is the noise in the pixel; and Pifip 
gives the total contribution from the foreground sources under consideration. We take the CMB signal and 
the foreground template to be already smeared by the beam and any other instrumental effects: 

s p = J d 2 x B(x,x p )T(x) 
= J d 2 x B(x ■ Xp)T(x) 

= ^2B e a hn Y em (x p ) , (2) 

tm 

where B(x, x p ) gives the response of the beam at position x when pointed at pixel p. In the second equality 
we assume that the beam is azimuthally symmetric around pixel p, and in the third equality we transform 
to spherical harmonics with indices (£, m). Here, Bt is the transform of the symmetric B(x ■ x p ) and ae m is 
the transform of T(x). For more information on working with asymmetric beams and pixels, see Wu et al. 
(2001); Souradeep & Ratra (2001). 

Note that this model, Eq. 1, is essentially the same as that considered in Stompor et al. (2002), although 
that work considers a broader class of templates, fi P , allowing them to be associated with an amplitude 
synchronous with instrumental characteristics rather than the sky. That work then derives techniques, 
analogous to the ones described below, to determine or marginalize over these amplitudes. 

Under this model, we can ask two separate questions: 

1. What are the foreground amplitudes, (5{l 

2. What is the underlying signal, s p ? Perhaps more important is the related question: What is Cp. That 
is, what are the statistics of the cosmological component, s p , taking into account the possible presence 
of foreground contamination? 
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To answer these questions we must first assign appropriate distributions to the quantities in Eq. 1. As 
usual, the noise is taken to be distributed as a zero- mean Gaussian with covariance N pp > = (n p n p i) (assumed 
known beforehand, although we can also apply iterative techniques (Ferreira & Jaffe 2000; Dore et al. 2001) 
to determine it simultaneously with the CMB and foreground signals). This gives a likelihood function 



P(d\/3,s) 



1 



\2irN\ 1/2 
1 



exp 



(d-s-pffN'Hd-s-pf) 



(3) 



where we have left off indices and used matrix notation. 



In order to determine the underlying signal power spectrum, we will need to assign an appropriate 
Gaussian distribution, with variance given by 



jpp' 



(s P s p >) = - BtC e P e (cos6 pp <), 



(4) 



where C\ = (\ag m \ 2 ) is the underlying CMB power spectrum, 9 pp > is the angle between pixels p and p' , and 
the Pi are the Legendre polynomials. This gives a distribution 



P(s|c,) = ^ exp 



(5) 



All of this is as in the usual foreground-free case. The likelihood for the data d (given (3 and s) is a 
Gaussian, now with mean '^2 i Pifi P + s p . Finally, however, we must assign a prior to the amplitudes 
To remain sufficiently general, we will allow an arbitrary Gaussian for each amplitude, with 1 (/?,) = 0, and 
(/3 2 ) = a%. The prior is then 



p<jt) = U- 



■. exp 



i A 2 



2^ 



(6) 



We can take erg, — > oo to give a 'non-informative' distribution. This is algebraically easier to deal with than 
the equivalent unbounded uniform distribution. 

If we combine these priors with the likelihood, we use Bayes' theorem to get the posterior distribution 



P{Pi,8 p \C e , d p ) oc P{p t )P(s p \Ci)P{d p \p t , s p , d) 



(7) 



where P((3i) and P(s p \Ci) are priors, and P(d p \/3i,s p ) is the likelihood. Then, in order to answer each of 
the above questions, we marginalize over [3 to get the posterior for s; s to get the posterior for /?; and finally 
marginalize over both, after giving s the appropriate prior variance appropriate for a given C(. Because of the 
linear, Gaussian form of our likelihood and priors, all of these marginalizations essentially involve "completing 
the square," and each of the posteriors remains a Gaussian distribution. For the same reason, these results 
can also be derived on the usual minimum-variance grounds rather than in this Bayesian formalism. 



: If desired, we can include a known mean for the foreground contribution by simply subtracting it from the data at the 
outset. 



- 5 - 



First, we calculate the posterior distribution of (3 by marginalizing over s p : 



P(P\d P ,C e ) 



ds P(0,a\d,C t ) 



1/2 



exp 



This is a Gaussian with mean 



and covariancc 



A = (Pi) = [f{S + N)- 1 /] J (f T (S + Ny'd), 
= ((/3 - 0i)(/3 - h)) = [f(S + N)- 1 /]^ 



(8) 



(9) 



(10) 



We have already taken the prior variances a pi — > oo in this expression. As is often the case in linear problems 
with normal errors as we have here, this is just the least squares solution for the amplitudes, (3{. 



The posterior for the signal, s p , is also a Gaussian. It has mean 

s p = S(S + N + a 2 pff T )~ 1 d 

and covariance 

({s-s p )(s-s p ,)) = S(S + N + alff T )-\N + alff T ) 
= S-S(S + N + a a p ff T )- 1 S. 



(11) 



(12) 



This gives us the Wiener Filter as the mean of this distribution. Note that we have left in a finite prior 
variance for We can take these prior variances to infinity using the Sherman-Morrison- Woodbury formula, 
which states 

(W + JBfy 1 = W- 1 - W-V(/ T W -1 / - B- 1 )- 1 / 7 ^- 1 . (13) 

Setting W = S + N and B^ 1 = oj^l — > here is equivalent to marginalizing over our template amplitudes. 
Note that, as in Stompor et al. (2002), this same formula can be applied in the determination of the original 
map, d p , in order to marginalize over such modes at an earlier stage. 

Finally, we wish to determine the power spectrum, Cg. We do this by starting with the expression in 
Eq. 7, and using Bayes' theorem yet again: 



P(C e \d p ) oc P(Ci)P(d p \C t ) - P(C t ) J dp dsP(/3, s\d, C t ) . 



(14) 



(Equivalently, we could have considered Cg a "parameter" from the start, and just given it a delta function 
prior when calculating the distributions of (3 and s.) In this case, it is worth writing out the entire expression: 



P(C e \d) = P(C e )- 



■ exp 



2 



|2ttM| 1/2 

where the covariance matrix is given by 

M pp , = (d p dp,) = ^ fi P (0if3i')fi'p' + (spSp') + (n p n pl ) 



(15) 



— ( a pf T f)pp' + Spp' + n pp' 



(16) 
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Recall that the signal covariance is a function of the underlying power spectrum, Ci, as in Eq. 4. We then 
maximize this with respect to the Ct (or bands thereof with known shape) as in Bond et al. (1998) as 
implemented in MADCAP 2 (Borrill 1999). In this expression, we have kept a finite variance for the dust 
prior, cr 2 ,. In the absence of a known dust contaminant spectrum (or to be most conservative), we can 
take up — > oo. The numerical implementation of this can be done using the Sherman-Morrison- Woodbury 
formula, implemented in the current version of MADCAP. 

Note that it is the linear nature of our model Eq. 1, along with our assignment of Gaussian priors and 
likelihoods, that allows these analytic simplifications. This makes it somewhat more difficult to determine 
or marginalize over more physical parameters such as, say, the spectral index of the foreground spectra. 
However, we can use this formalism to determine several /3,, each corresponding to a different region of the 
sky — allowing us to take account at some level of spatial variations in the spectrum. For the present case, 
however, the final signal-to-noise is too low for this to be fruitful. 

In the Bayesian picture, our method simply ignores any information associated with a pattern on the 
sky matching our chosen templates. Thus, even if our template were "wrong" (i.e., if they did not accurately 
reflect the pattern of foreground emission on the sky), the final estimate of the power spectrum, including 
its error bars, would be no less correct than that without applying this method at all. In the simplest case, 
the maximum likelihood power could decrease, but the error bars would increase to take this into account. 



2.1. Noisy Foreground Templates 

So far, we have only considered foreground templates, /, that are accurately known. In many cases, 
however, the pattern will also be the result of a noisy measurement. That is, we supplement our model of 
the data, Eq. 1 with a new set of pixel measurements, 



(17) 



where e p is the new data, v p is noise, satisfying (v p v p i) = E pp >, and we can again assign a Gaussian error 
distribution for v p . Now, to proceed we must multiply this new distribution by the likelihood of Eq. 3 or the 
full Eq. 7, and marginalize over the now unknown template, /: 



P(de\/3,s) 



\2irN + [3 2 E\ 1/2 
1 



exp 



(d-s- (3e) T {N + f3 2 E)- 1 {d - a - pe) 



(18) 



After marginalization, the variance is increased by /3 2 E over the case where the foreground template is known 
exactly. As before, we can also put in a prior for the signal amplitude, s (Eq. 5), and marginalize, giving 
the equivalent of Eq. 15, 



P(C e \de) = P{C t ) 



\2n(S + N + (3 2 E)\ 



1/2 



exp 



-(d - pe) T (S + N + f3 2 E)- 1 (d - j3e) 



(19) 



In both cases, the marginalization over the foreground template has destroyed the linear nature of the 
problem, and the unknown amplitude, (3 now appears in the effective correlation matrix; we can no longer 
find analytic formulae to determine nor marginalize over the amplitude 



2 http://www.nersc.gov/~borrill/cmb/madcap. html 
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In fact, the problem for a noisy template is formally identical to the cross-calibration of cross-comparison 
problem, as discussed in Knox et al. (1998): the amplitude, j3 is just the relative calibration of the foreground 
contribution to the two datasets. As in that work, however, we must use somewhat more laborious numerical 
techniques to deal with the more complicated problem. 

3. Application: Diffuse foregrounds and MAXIMA-1 

Here we apply these techniques to the data from the MAXIMA- 1 experiment; the hardware is described 
in Lee et al. (1999) the resulting maps and power spectra are described in Hanany et al. (2000); Lee et al. 
(2001), the cosmological implications are presented in Balbi et al. (2000); Stompor et al. (2001); Jaffe et al. 
(2001) and details of the data analysis in Stompor et al. (2002). In this work, we consider maps made from 
the data of individual photometers at various frequencies: three at 150 GHz, two at 240 GHz, and two at 
410 GHz. Note that only a subset of these maps were used for the cosmological results presented in the 
previous MAXIMA-I papers. A full account of the systematic checks on the MAXIMA-I data and on the 
consistency between channels is given in Stompor et al. (2003). 

The primary aim of the MAXIMA experiment is to observe the cosmological component of the CMB. 
The MAXIMA-I field was thus explicitly chosen to be low in Galactic foreground emission. Nonetheless, our 
methods let us take advantage of the large number of pixels simultaneously, even though the dust contribution 
is insignificant in a single one. 



3.1. Dust Emission from Galactic Cirrus 

As foreground templates, we primarily work with the combined IRAS/DIRBE infrared dust maps pre- 
sented in Schlegel et al. (1998). These are extrapolated from 100-240^m to our CMB observation frequencies 
using dust emissivity models constrained by COBE/FIRAS spectral observations (Finkbeiner et al. 1999). 
One can use these emissivity models to compute a dust temperature and column density in each pixel from 
the DIRBE 100 and 240^m maps, and extrapolate to much lower frequencies. Although our formalism could 
be used to extrapolate from the observed infrared frequencies down to our observations at 100-400 GHz, 
neglect of the dust temperature variation from pixel to pixel can cause errors of a factor of ~ 2, so we use 
the FDS99 models that already take account of this variation explicitly. Eight models are computed; four 
single-model components (power laws) and four two-component models (where the two components have a 
fixed mass ratio, and have temperature coupled to the radiation field in a self-consistent way). The param- 
eters of the model are given in Table 1, adapted from FDS99. The 4 parameters describing each model are: 
a\, a2, the emissivity power-law indices for components one and two, /i, the fraction of power absorbed and 
re-emitted by component 1 (unrelated to the foreground template vector / above), and q\/q2 where qi is the 
IR/optical opacity ratio for component i. The two-component models listed in the table are: 5) T\ = T 2 , 
indices from Pollack et al. (1994), 6) ct\ = a 2 = 2 as in Reach et al. (1995), 7) like 5 except temperatures 
can float, and 8) all four parameters floating. Because 8 gives the best chi-squared, it is the preferred model, 
although 7) is not appreciably different. Here, we will concentrate on their overall best-fit model 8, with 
on = 1.67, a 2 = 2.70, /i = 0.0363, and q 1 /q 2 = 13.0. 

In Figures 1,2, and 3 we show the MAXIMA-1 data as well as the SFD extrapolated dust maps (using 
their preferred Model 8) at each of these frequencies. Note that the temperature scale for the dust at 150 and 
240 GHz is stretched considerably compared to the data - the expected RMS dust contribution is ~ lfiK at 
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150 GHz, compared to the ~ 3C%K (signal plus noise) RMS of the 150 GHz CMB map. 





Fig. 1. — Left: MAXIMA-I map of microwave emission at 150 Ghz. Right: Map of dust emission, extrapo- 
lated from IRAS/DIRBE maps to 150 GHz. Colorbars give the temperature in /j,K. Note that the difference 
in the temperature range between the total and dust emission is a factor of ~ 50. 



Because of this large difference in the amplitude of CMB and dust emission, we do not expect to be able 
to see the dust contribution to a single pixel in the maps. However, our procedure for estimating (3 gives 
us the usual (A^dof) 1 ^ 2 advantage when considering the whole map (Adof gives the number of degrees of 
freedom in the map, equal to the number of pixels less any degrees of freedom marginalized over in making 
the map or by the methods described here). We can further increase the signal-to- noise by combining the 
individual photometers at a given frequency. This allows us to compare the different models offered by SFD, 
as detailed in their work. 

In Figure 4, we show the observed amplitude (3 for each SFD model, for the MAXIMA-I detectors 
combined at each of 150, 240 and 410 GHz, and combined over all frequencies (the latter makes sense only 
if the overall spectral shape of the model is correct over this frequency range). The data do not strongly 
prefer any single model. However, a few results are evident. 

The dust signal is not strongly detected at 150 or 410 GHz for any model. That is, at these frequencies 
an amplitude of /3 = is not disfavored. Indeed, some of the models are disfavored at roughly the one sigma 
level from the 240 GHz data, with the one-component model 1 the most disfavored. At 410 GHz, the dust 
signal is stronger, and (3 = 1 (the model prediction) is preferred over (3 = in all cases. 

Despite the only marginal preference for a non-zero signal, these results do contain important informa- 
tion. Even in the cases where (3 = is allowed, the results can be interpreted as an upper limit on the dust 
amplitude at these frequencies and in this area of sky. As is evident from equations 9 and 10, the observed 
amplitude and error scale inversely with the template amplitude. Thus, although (3 = 1 is acceptable for all 
of these models, models predicting dust emission a factor of a few higher at any of these frequencies would 
be strongly disfavored. 

In Figure 5, we show the results for SFD's Model 8 — their overall best fit — in detail. We show each 
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Fig. 2. — Left: MAXIMA-I map of CMB Emission at 240 Ghz. Right: Map of dust emission, extrapolated 
from IRAS/DIRBE maps to 240 GHz. Note that the difference in the temperature range between the total 
and dust emission is a factor of ~ 17. 




14. S 15.0 15.2 15.4 15.6 15. B 1S.0 14. B 15.0 15.2 15.4 15.6 1 5.E 16.0 

RA (hours) RA (hours) 



Fig. 3. — Left: MAXIMA-I map of CMB Emission at 410 Ghz. Right: Map of dust emission, extrapolated 
from IRAS/DIRBE maps to 410 GHz. Note that the difference in the temperature range between the total 
and dust emission is a factor of <~ 40. 
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Fig. 4. — Each panel gives the observed amplitude of the dust template for the model number given on the 
horizontal axis at a given frequency. The top panel averages over all frequencies. 
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of the individual detector amplitudes, the frequency averages, as well as an overall average. In Figure 6, 
we rescale the observed coefficients (i to give the actual values of the observed dust emission amplitude in 
thermodynamic temperature units. 

This figure emphasizes that the techniques described here are not only useful for the removal of CMB 
foregrounds, but are more generally useful for estimating the foreground spectrum and extrapolating it to 
regimes where it may be completely negligible in an individual pixel. 

3.2. Synchrotron Emission 

We can of course use the same algorithm with other sources of foreground emission for which we have 
maps. The FSD code also provides an extrapolation synchrotron emission as measured by Haslam et al. 
(1981); Reich & Reich (1986); Jonas ct al. (1998). 

These surveys were reprocessed, destriped, and point-source subtracted (D.P. Finkbeiner 2002, private 
communication) and are available to the public as part of the dust map distribution 3 . 

The Haslam survey is full-sky, the Reich & Reich survey is in the north, and the Rhodes survey in the 
south, so for every point on the sky at least two frequencies are available. The surveys are beam-matched (at 
a one degree resolution) and used to determine a power law for each pixel on the sky. This power law is then 
extrapolated to the frequency of observation. Because the synchrotron spectrum is known to fall faster than 
a power law at high frequencies, this prediction should be interpreted as an upper limit on the synchrotron 
emission. If the 408 - 2326 MHz maps are contaminated by significant free-free emission, the power law 
slope is shallower than it should be, making it even more of an upper limit. The fact that even this upper 
limit is smaller than 3/xK at 150 GHz over the MAXIMA area implies that the synchrotron emission should 
be totally undetectable. Figure 7 confirms this to the extent possible with this data. 

3.3. The CMB power spectrum 

Now that we have measured the overall level of dust emission in the MAXIMA-I field, we can now ask 
the other questions posed in Section 2: what does the CMB itself look like? In figure 8, we show the power 
spectrum of CMB temperature fluctuations from the combined 150 GHz photometers (as used in Lee et al. 
(2001) and Stompor et al. (2001)). One set of points shows the spectrum ignoring the contamination from 
dust and synchrotron emission (073 — > in equations 15 and 16); the other set marginalizes over the dust 
emission with the known morphology of FSD's model 8 (erg — > 00). We see that the marginalization has 
very little effect — much smaller than the error bars. This is consistent with our knowledge of the power 
spectrum of high-latitude dust emission, with Ce ~ £~ 3 (or perhaps closer to l~ 2 in some parts of the sky), 
dominating only at the largest scales. This implies that the spatial pattern of the dust (at least in the 
MAXIMA-I patch) is essentially incompatible with that of an isotropic Gaussian field on the sky. 



3 http://astro.berkeley.edu/dust 
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Fig. 5. — Ratio of expected amplitude of dust emission (SFD Model 8) to that observed in the MAXIMA-I 
maps, as a function of detector frequency (offset slightly for legibility). Thin error bars with circles are 
individual detectors, thick error bars with crosses at the individual frequencies (150, 240, 410) combine 
these, and the single error bar at / ~ 300 gives the overall average. 
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Fig. 6. — Thermodynamic temperature of dust emission observed in each of the three MAXIMA- 1 frequencies. 
The curve gives the average emission from Finkbeiner et al's Model 8 predicted in the MAXIMA-I observation 
area. 
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Fig. 7. — Observed synchrotron emission relative to that expected, as in figure 5. 
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Fig. 8. — CMB Power spectrum ignoring the effect of dust and synchrotron contamination (left crosses) and 
marginalizing over it (right circles). 
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4. Discussion and future applications 

We have derived a technique for measuring and accounting for the effect of foreground emission on CMB 
observations, for the case where the morphology of the contaminant is known, but when its spectrum is un- 
known or imprecisely measured. We have applied these techniques to the MAXIMA-1 data and observations 
of dust and synchrotron emission. The dust emission in the MAXIMA-I region is consistent with models and 
observations at higher frequencies, although it has negligible effect on the measured CMB power spectrum. 

The MAXIMA-1 observing region was specifically chosen to be a region of low dust contrast. As high- 
resolution CMB observations cover more of the sky with higher sensitivity, these techniques will become 
more important for the separation of the various components. 

These techniques have many further applications, some alluded to above. We can allow for spatial varia- 
tion in the foreground spectrum and/or inaccuracies in our foreground templates. Most straightforwardly, we 
can apply the technique separately to individual patches (with, say, different dust temperatures) and allow 
the foreground amplitudes to float separately between them. Such a technique could be applied iteratively 
on smaller patches until the signal-to-noise of the result decreased too far. In particular, we would certainly 
split the full sky up into regions within and outside of the Galactic plane, where we know the foreground 
properties to differ. Other prior knowledge will affect the algorithm in different ways. If we thought that 
the foreground spectrum was approximately a power law, f p oc is b , we can estimate (or marginalize over) the 
power law index. The change in the foreground amplitude is fpdblnv, equivalent to the (3f p in Eq. 1. 

More ambitiously, if we have knowledge of foreground polarization, these techniques carry forward 
identically, although such measurements may not be readily forthcoming. 

Finally, we have discussed here the case where the foreground template is known with considerably more 
accuracy than the CMB measurement. As CMB observations are performed with higher sensitivity, we will 
need to deal with foreground templates with errors of their own. Although the calculations and resulting 
algorithms are considerably more complicated, the same general setting can be used; the results are similar 
to the cases discussed in Knox et al. (1998). 
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